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Abstract 

In close analogy to fermionic many-body theory the truncation of the bosonic 
BBGKY density matrix hierarchy on the two-body level leads to a coupled set of 
nonlinear equations of motion for the one-body density matrix and the two-body 
correlation function. These equations provide a nonperturbative description of the 
nonequilibrium time evolution of particle number conserving bosonic many-body 
systems including the dynamical resummation of parquet-like diagrams. Within this 
framework we study the properties of a pionic system as a function of temperature 
and density with focus on two-body quantities. For each temperature we find a 
related pion density for which the relative strength of the two-body correlation 
function assumes a maximum and the pionic system is far from the mean-field 
limit. Since these correlated phases up to T=200 MeV only appear at rather low 
pion density, the hot and dense pion gas as generated in ultrarelativistic nucleus- 
nucleus collisions should be well described within mean-field theory; i.e. the HBT 
analysis of pion sources from vr — vr correlations should remain valid even in the case 
of strongly interacting pions. 
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1 Introduction 



The theoretical description of strongly interacting hadronic many-body systems in general 
requires nonperturbative methods. A systematic way for constructing a wide class of such 
nonperturbative approaches is given by the truncation of Green function hierarchies as 
obtained by inserting the cluster expansion in terms of the corresponding connected Green 
functions. In the lowest order truncation scheme one then ends up with the mean-field 
level, i.e. TDHFQ The next order truncation scheme, i.e. neglecting the connected six- 
point function (the three-body correlation function) as well as all higher order correlation 
functions leads to the equations of two-body correlation dynamics, which in the fermionic 
case have been denoted by NQCD^ The NQCD-method has already been successfully ap- 
plied to the most important fermionic many-body problem in hadronic physics, the nuclear 
many-body problem It has been shown to guarantee a simultaneous resummation 

of ring- and ladder-diagrams in vertical and horizontal direction (a parquet resummation) 
and thus to adequately take into account long- as well as short-range nucleon-nucleon 
correlations p. 

It is the purpose of the present paper, to develop a related nonperturbative theory 
for bosonic systems. The most important bosonic many-body problem on the hadron 
level is that of a pion gas, which can be generated e.g. in an ultrarelativistic heavy- 



ion collision ||T0|-||T2|. Here especially the Hanbury Brown-Twiss (HBT) analysis of pion 
sources from tt — vr correlations |T3[-|2^ might become questionable due to the strong 
interaction between the pions. In this work we will apply our method to study the 
properties of two-body quantities in such systems. 

The paper is organized as follows: In Sect. ^ we will derive the equations of motion of 
two-body boson correlation dynamics starting from the Heisenberg equations of motion 
for an arbitrary, particle number conserving many-body Hamiltonian with a two-body 
interaction. In Sect. || we will investigate the topological structure of the equations and 
analyze the various interaction terms diagrammatically. 



Sect. 4.1 to Sect. 4.4 are devoted to a presentation of numerical results for pionic 



systems, i.e. the comparison of different limiting cases for the resummation of diagrams 
(Sect. [4.1| ), the computation of correlated near-equilibrium states (Sect. [4.2| ), the effect 
of correlations on two-body quantities (Sect. |i7^ ) and signatures for a phase transition 
of the system as a function of density and temperature (Sect. [4.4|) . We will shift the 
specification of the Hamiltonian density for interacting pions as well as the final explicit 
equations of motion to the Appendices and 0, respectively. 



2 The equations of motion 

In this Sect, we develop the mathematical apparatus of correlation dynamics for bosonic 
many-particle systems described by an arbitrary, particle number conserving Hamiltonian 
with a two-body interaction. 

In this respect, we first introduce the formalism of density matrices: let and 
be the creation and annihilation operators for a system of bosons, where a stands for a 
complete set of quantum numbers characterizing an element of an arbitrary, orthonormal 

^Time Dependent Hartree Fock 
^Nuclear Quantum Correlation Dynamics 
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basis of the corresponding one-body Hilbert space, with canonical equal-time commutation 
relations 

K, = 5af3 , [a^, a/s] = [a^, 0)3] = . (2.1) 

The matrix elements of the n-body density matrix are then defined by 

(Pn)ai...a„a;...< = (o^, • • .olj^ a„„ . . .a„, ) , (2.2) 

i.e., by the expectation value^ of the normal-ordered operator product of the correspond- 
ing creation and annihilation operators, all of these considered at equal-time t in the 
Heisenberg picture. The matrix elements of the density matrices can thus be considered 
as equal-time Green functions. 



From the canonical commutation relations (|2.1|) we immediately obtain the symmetry 
relations 

(Pl)aQ' = (Pl)^'^ 5 (P'i) af3a' P' ^ (P'^) a' iS'aP ' 

iP^)af3a'l3' = (P^) f3aa' f3' = iP'^) aPP' a' = (P2)/3a/3'a' • (2-3) 

Analogous relations hold for all higher density matrices. 
We now consider a Hamiltonian of the general form 

with a pure two-body interaction v. From the hermeticity of the Hamiltonian we get 

taa' = Oa , {a(3\v\a' (3') = {a'(3'\v\a(3)* . (2.5) 

The time evolution of an arbitrary operator O with no explicit time dependence is given 
by the Heisenberg equation 

idtO = [O, H] , (2.6) 

where dt denotes the total time derivative. Then the time derivative of the operator 
product 

n n 

(^a'---^a'^a„---(lai = «q,' • • -^q,' «ai • • = II '^-q' H O-ofc , (2.7) 

i=l k=l 

which defines the n-body density matrix (|2.2|) via the expectation value, reads: 



idt 



n< n 

L \i=l / \k=l 



E E ( n V ) I n 1 «A - txa'^ I n v. 1 4 ( n « 



j=l A I \i=l ■/ \k^j=l I \i^j=\ ^ I \k=l 



Oik 



+ ^ E E E {(^j(^k\v\XiX2)s n 4; n 1 ^AiflAa 

^j=lfc=lAiA2 \i=l / \ljtj,k=l 



•^in general an ensemble average with an arbitrary nonequilibrium quantum statistical density operator 
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n 



i=lfc=i+lAiA2 \ii^i,k=l 7 \l=l 



1 ^ 



n \ n 



^j=lAiA2A3 [ \i=l / \Mi=l 



n 

-(AiAalt^lAgaj)^ ( H 4j ^i^L (^H ^^^j 1 (2.8) 
with the symmetrized two-body matrix elements 

{aj\i\v\\2\'i)s = {aj\i\v\\2\3) + {aj\i\v\\s\2) ■ (2.9) 

In the middle two terms one can use either Ylk=i J2k=j+iJ because the expressions in the 
sums are symmetric in i and j. The equations of the BBGKY density matrix hierarchy 
now follow by taking the expectation value on both sides of (p.8|)p|. With ( p.2| ) we obtain: 

n 

'^dt{pn)ai...ana[...a'„ — E] 5Z {tajX{Pn)ai...aj^i\aj+i...ana[...a'„ 
i=l A 

~^>^oi'. {Pn)ai...a„a[...a'._j^\a'.^j^...a'^^ 
n j-l 

)Q:i...afc_lAiQ:fc+i...Oj_iA2«j + i...«ri«'i-..aJ, 



1 ^"-^ r 
+ 1^YYY1 {{ajak\v\XlX2)s{Pn)c 



2 j=l k=l A1A2 

~(^1^2|^|ajafc)5(p„)Qi...a„a;...a;^_^Ai«'^_^^...a^._jA2a;+i...<} 
1 " 

+ X! {{(^jXl\v\X2X3)siPn+l)ai...aj.iX2X3aj+i...a„a[...a'^\i 

^ i=l A1A2A3 

-(^l-^2|^^|A3aj)5(p„+l)ai...a„A3a;...a^_iAiA2a;+i...<} • (2-10) 

The time evolution of the density matrices consequently is given by a coupled system of 
equations of first order in time, where the equation of motion for p„ couples to p„+i. For 
practical purposes the hierarchy ( |2.1(j| ) has to be truncated. 

For the special cases of pi and p2 - which we will examine furtheron - the equations 
read explicitly: 

^dt{p\)aoc< = Y {^»xipl)xa' - txa'ipl)ax} 
X 

+ ^ {(a^l|^l^2A3)5(p2)A2A3a'Ai " (Al A2 1^; |a'A3)5(p2)aA3AiA2 } (2.11) 

A1A2A3 

and 

idt{p2)af3a'l3' = Y ax{P2) X/Ba' P' + i f3x{.P2) aXa' P' — i Xa' {pi) a^XP' — txp' {.P2) ajia' x} 
X 

+ \ Y {(a/^l^l^l^2)5(p2)AiA2a'/3' - (Al A2 |a'/3%(p2)a/3Ai A2 } 
A1A2 

+ ^ Y {("-^l|^|-^2A3)5(p3)/3A2A3a'/3'Ai + (/^Ai A2A3)5(p3)„A2A3a'/3'Ai 
A1A2A3 

-(AiA2|t;|A3a%(p3)a/3A3/3'AiA2 " ( Al A2 |f | A3/5%(p3)a/3A3a'Ai A2 } • (2-12) 



*One can therefore regard the BBGKY hierarchy as an operator identity. 
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A suitable truncation scheme for (|2.1CI|) is based on the cluster decomposition of Green 
functions, i.e. on the fact that any n-point Green function can be decomposed into a sum 
of products of connected Green functions of equal or lower order. The truncation then 
is performed by neglecting all connected Green functions higher than a certain order, i.e. 
in our case all connected Green functions of higher order than the four-point function (or 
two-body correlation function) This strategy is based on the assumption that the 

connected (correlated) parts of the corresponding Green functions become increasingly 
unimportant at higher order |P, p3| . 

A reduction to the two-body level requires the equations of motion for pi ( |2.11|) and p2 



( p.l2| ), which contain all density matrices up to p^. Thus we need the cluster expansions 
of pi, p2 and ps up to the two-body correlation function. With {pi)aa' = Paa' and c„ 
denoting the n-body correlation function we obtain: 

(Pl)ao' = Paa' , (2.13) 

{P2)al3a'f3' = Paa'PplS' + Pa/S'PfBa' + {c2)aPa'(i' 

= (p2o)o/3o'/3' + (c2)a/3a'/3' , (2-14) 

(P3)a/37a'/3'7' = Paa' PlS/B' P-yY + Paa' P/Sy P-yl3' + Pa/B' P/Sa' P-yY 
+Pa/3'Pl3YP-ya' + Pa-/' P^/3' P-ya' + Pa-y' P^a' P-y/S' 
+Paa'{c2)f3y/3'-y' + Pa/3' {C2) f3^a'Y + Pay' {C2) i3-y/3' a' 
+P/3l3'{c2)a-ya'-y' + Pf3a' {C2) ay/B'Y + P/Sy {C2) a^a' 13' 
+P^'y'{c2)af3a'f3' + P^a' {C2) 0/3/3' Y + P-y/3' {C2) a/3a'y 

a/3^a' /B'-y' 
= {P3o)a/3^a'/3'Y 

+ (1 + "Pa/B + 'Pa'y){^ + 'Pa'/3' + 'Pa'y) Paa' {C2) /3-y/3'y 
a/B-ya' /3''y' 1 

(2.15) 

where Va/3 is the permutation operator interchanging the indices a and /?; 

{P2o) a/3a' /3' = (1 + 'Po/b) Paa' P /3/3' (2-16) 

and 

(P3o)Q:/37a'/3'7' = (1 + Pa/3 + Pa'y){^ + P/3-y) Paa' P/3/3' P^Y (2-17) 

are the uncorrelated parts of the two- and three-body density matrices. 

The explicit expressions for the cluster expansions can be derived from the generating 



functionals of full and connected Green functions p4 . 

In our present case we make the additional assumption, that operator products with 
different numbers of creation and annihilation operators, i.e. operators not conserving 
particle number, have vanishing expectation values, e.g. 

{alai) = {a^ap) = . (2.18) 
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Neglecting C3 in (|2.15|) and inserting the truncated cluster expansions into the equa- 
tions of motion ( |2.11| ) and ( p.l2| ), we arrive at a closed system of coupled nonlinear dif- 



ferential equations for p and C2. In order to compactify these lengthy equations (in close 
analogy to fermionic correlation dynamics [0, we introduce the following abbreviations: 

Ugg' = PA1A2 (Q^A2|f |a^Ai) (mean — field) (2.19) 

Ai A2 

haa' = taa' + Uaa' (mean — field — Hamiltonian) (2.20) 

Qada'p' = Saa'Spp' + Paa'Spp' + ppp'S^a' (vcrtical Bosc - factor) (2.21) 

Qipa'p' = Ppa'^ais' " PafS'S^a' (horizontal Bose - factor) . (2.22) 

Since no higher correlation functions are included, C2 is simply denoted by c furtheron. 
The compactified equations of motion for p and c then read: 

idtpaa' = Y i^aXPXa' - hxa'Pax) (2.23a) 
A 

+ J2 (("Ai|?7|A2A3)cA2A3a'Ai - (AlA2|t;|a'A3)Cc,A3AiA2) , (2.23b) 
A1A2A3 

and 

idfCaiSa' 13' = Y i^aXCx/Sa' 13' + hpxCaXa'fi' " hxa'CafiXfi' — h Xl3' Ca/Sa' x) (2.24a) 
A 

+ Y {^3Xa\v\XiX2) [Q~i3X3X4iP20)xiX2a'f3' - QxiX2a'f3'iP20)apX3X4} (2.24b) 
A1A2A3A4 

+ J2 {^3X^X1X2) {Q^pXsXiCx,X2a' (3' -QZx2a'l3'Caf3X;X^} (2.24c) 
A1A2A3A4 

+ (1 + V^^) (1 + V^>^-) J2 {>^3>^M>^l>^2)sQL,a'X,CX2(3X,(3' • (2.24d) 

A1A2A3A4 

The approximation of the bosonic BBGKY density matrix hierarchy given by (|2.23| ) and 
( p.24| ) will be denoted as NCBCD (Number-Conserving Boson Correlation Dynamics) 



furtheron. The equations conserve e.g. particle number (A^) = J2a Paa and total energy 
(if) for a Hamiltonian with no explicit time dependence. 

We note in passing that the compactified bosonic NCBCD equations are formally very 
similar to the compactified NQCD equations for the nonrelativistic description of fermionic 
nuclear systems 0] with Bose-factors instead of Pauli-blocking-factors and symmetrized 
instead of antisymmetrized matrix elements of the two-body interaction and the two-body 
density matrix. 

The NCBCD equations guarantee a dynamical, nonperturbative resummation of ring- 
and ladder- diagrams in both vertical and horizontal direction and of the corresponding 
mixed diagrams without double-counting. Topologically this corresponds to considering 
parquet-like diagrams. 
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3 The diagrammatical structure of the equations of 
motion 

In this Sect, we investigate the topological structure of the NCBCD equations of motion 
for the pionic model specified in Appendix ^ and give a graphical representation of the 
terms appearing in the equations in order to illustrate the classes of diagrams that are 
resummed. 

The interaction v (see Appendix ^) is split up in the t-, u- and s-channel part: 

{ap\v\a'p') = {a(3\vt^^\a'(3') + {a(3\vs\a'(3') . (3.25) 

In order to allow for a unique identification, we now label the individual terms in the 
equations of motion: 

idtPaa' = (tciXPXa' - txa'Pax) (3.26a) 
A 

+ J2 (UaXpXa' - Uxa'Pax) (3.26b) 
A 

+ Y (("-^i|^|-^2A3)ca2A3q'Ai - (AlA2|t;|«'A3)CaA3AiA2) , (3.26c) 
A1A2A3 

idtCa(5a'P' = X! i^aXCxpa'P' + tf^x^aXa'P' — ^Xa'CapXp' — t X/3' Cafia' x) (3.27a) 
A 

+ iUaXCXfia'P' + UpxCaXa'l3' — Uxa'Caf3Xf3' — Uxp' Cafia' x) (3.27b) 
A 

+ Y (A3A4|^'t,,i|AiA2) (Q~^A3A4(P2o)AiA2a'/3' - <5AiA2Q'/3'(P2o)a/3A3A4) (3.27c) 
A1A2A3A4 

+ Y (A3A4|fs|AiA2) (Q~^A3A4(P2o)AiA2a'/3' - QAiA2«'/3'(P2o)a/3A3A4) (3.27d) 
A1A2A3A4 

+ Y (A3A4|^^t,«|AiA2 )(Q CXiX2a'/3' - QxiX2a'l3'(^al3X3X^ (3.27e) 

A1A2A3A4 

+ Y (A3A4|^^s|AiA2) (Q~^A3A4CAiA2a'/3' - QAiA2Q'/3'Ca/3A3A4) (3.27f) 
A1A2A3A4 

+ {l+V^p){l+V^>p>) Y {h>^M,u\>^l>^2)sQix,a'X,CX2PX,(3' (3.27g) 

A1A2A3A4 

+ (1 + P^^) (1 + P,,^0 J2 (A3A4|^s|AiA2)5QiA,a'A3CA2/3A4/3' ■ (3.27h) 

Ai A2A3A4 



3.1 The one-body equation 

The term (3.26a) accounts for the free propagation of the one-body density matrix. In 
general all terms are of hermitean structure, i.e. for each process there is the conjugate 
process with the incoming and outgoing indices interchanged. 
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The term ( p.26b ) accounts for the influence of the unrenormahzed mean-field on the 



propagation of the one-body density matrix, (3.26a) and ( p.26b| ) together then are equiv- 
alent to the TDHF approximation. 

In order to investigate the structure of ( p.26b| ), we first examine the part with positive 
overall sign: 

A 

E PAiA2("A2|t;t,«|A3Ai)pA3a' (3.28a) 
Ai A2A3 

+ E PAiA2(ttA2|t;t,«|AiA3)pA3a' (3.28b) 

Ai A2A3 

+ J2 PAiA2(ttA2|Ws|A3Ai)pA3a' (3.28c) 
Ai A2A3 

+ E PAiA2("A2|t;s|AiA3)pA3a' . (3.28d) 
A1A2A3 

We now use the following convention: In all diagrams of this Sect, solid lines stand for the 
link between two equal indices or for external indices. They describe the propagation of 
the pions participating in the process. In the graphical representation of terms containing 
a cr-propagator we use a horizontal dashed line for a spacelike and a vertical dashed line for 
a timelike a. At this point it is useful to note, that in momentum space Vt^u is proportional 
to a cr-propagator depending on the difference of two on-shell four momenta while Vs is 
proportional to a a-propagator containing the sum of two on-shell four momenta, where 
for Vs (up to the formfactor) the factor multiplying the propagator is only half of that for 
^i,tJ3 (see Appendix 

The contributions (3.28a) and (|3.28c|) are depicted in fig. ^ and correspond to the 
Hartree approximation. The terms ( |3.28b| ) and (|3.28d| ) are depicted in fig. ^ and corre- 
spond to the exchange terms in the Hartree-Fock approximation. The graphical repre- 
sentation of the conjugate terms is obtained by interchanging the incoming and outgoing 
external indices. 

The coupling of the equation of motion for the one-body density matrix to the two- 
body correlation function is given by (|3.26cj ). In fig. |^ we show only one of the two 
conjugate processes. Both conjugate processes together have a gain-loss structure as it 
is known e.g. from collision terms in transport equations However, ( p.26c|) cannot 



directly be identified with a collision term, since it also renormalizes the mean-field. 

In order to present the last point more explicitly, it is useful to introduce a G-matrix 
and a selfenergy within the formalism of equal-time density matrices. The G-matrix 
Gai3a'f3' IS defined by 

E G^a/3AiA2(P2o)AiA2a'/3' = E ("/^l^l A1A2) (P2)aiA2Q:'/3' • (3.29) 
A1A2 A1A2 

Inserting ( ^.291 ) into the equation of motion for p leads to 

idtPaa' = E (taXPXa' " t\a' Pa\) + E (^f A3AiA2PAia'PA2A3 - G'f^;^2a'A3PQAiPA3A2) (3.30) 
A A1A2A3 

'"'The difference of two on-sliell four momenta of particles of equal rest mass is always spacelike, the 
sum is always timelike. 
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with G^f^^ip, = Gai3a'/3' + Gaf3/3'a'- The G-Hiatrix terms in ( p.30|) contain all corrections to 
the free propagation, i.e. also the renormalized mean-field. 

A further compactification of the one-body equation can be achieved by introducing 
the selfenergy according to 

^aa' = PAiAaG'^Aaa'Ai • (3.31) 

A1A2 



Inserting ( p.31|) into (|3.3CI|) gives 



idtPaa' — X! {i^aX + ^qa) P\a' — {txa' + ^Aa') Pax} ■ (3.32) 
A 

The renormalized mean-field can now be obtained as the hermitean part of T^aa', i-e. 

= 2 ^ (pAiA2G'f;,2a'Ai + PAiAa'^f'AaAi) = P^iAa i^^'^^) ^Xza' X^ ' (3.33) 

A1A2 A1A2 

i.e. as the selfenergy obtained from the hermitean part of the symmetrized G-matrix. 
The renormalized mean-field Hamiltonian can then be defined by 

^aa' = '^aa' + U^aai ■ (3.34) 

With the above definitions the equation of motion for the one-body density matrix reads 

idtpaa' - Y {KxPXa' " Ka'Pax} = {(ImJ:) ^XpXa' - paxilruj:) Xa'} , (3.35) 
A A 



where the term on the right hand side of (|3.35| ) corresponds to the collision term of a 



transport theory. The explicit derivation of a transport equation with a collision term 
containing the correct in-medium factors and transport coefficients from the NCBCD 
equations can be carried out in complete analogy to the fermionic case 0, ^ . 

3.2 The mean-field terms in the two-body equation 

The term (3.27a) in analogy to the corresponding one in the one-body equation accounts 
for the free propagation of the two-body correlation function. 

The mean-field part ( |3.27b| ) of the two-body equation has the same structure as in 



( p.26b|) ; for each of the external indices there is a diagram with a mean- field loop at the 
corresponding line. 

3.3 The Born-term 



The expressions (|3.27c|) and (|3.27d|) only contain matrix elements of the two-body in- 



teraction coupled to a Bose-factor and two one-body density matrices, but not the 
two-body correlation function itself; they describe exactly one elementary interaction pro- 
cess between two pions and thus correspond to the first Born approximation of scattering 
theory, i.e. ( p.27c| ) in the t- and u- channel and ( p.27d| ) in the s-channel. The two terms 
together are hence referred to as Born-term. 
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The scattering process is subject to an in-medium modification by the Bose-factor in 
the incoming or (for the conjugate process) in the outgoing channel. Using the single 
particle basis that diagonalizes p, i.e. 

Paa' = ^aa'na , (3.36) 

the summations over the two intermediate indices appearing in in ( |3.27d|) and (p.27c|) 
break down and we obtain factors oil + rta + np in the incoming channel and l + n'^ + n'^ in 
the outgoing channel; the Bose-factors thus cause an enhancement of the contribution of 
the corresponding process to the total amplitude proportional to the occupation numbers 
of the respective states. As an example, two of the diagrams generated by ( |3.27c[ ) and 
( p.27d|) are represented in fig. |^. 

3.4 Dynamical iteration in the vertical direction 

The terms (|3.27e|) and ( p.27j) each contain a matrix element of the two-body correlation 
function itself coupled to a matrix element of the two-body potential, leading to a dy- 
namical iteration of the interaction in "vertical" direction (i.e. along the time direction 
of the diagrams) . The same Bose-factor Q"" as for the Born-term appears in the incoming 
or the outgoing channel. The two terms with Q"" in the incoming channel are graphically 
represented in fig. |. 

The left diagram in fig. |^ shows, that the term ( p.27e ) leads to a dynamical resum- 



mation of ladder-diagrams in vertical direction, where in the intermediate states as well 
as in one of the external channels there is a Bose-factor . The right diagram in fig. 
^ shows, that the term (|3.27f|) leads to a dynamical resummation of ring-diagrams in 



vertical direction with the same factors of Q appearing as for the ladder diagrams. The 
first step of the resummation in each case is given by the Born-term; ( p.27e ) and ( p.27i|) 



together lead to a mutual dynamical iteration of the ring- and ladder-diagrams as shown 
in fig. |. 

In analogy to the fermionic NQCD case one obtains a Bethe-Goldstone equation for 
the G-matrix by neglecting ( p.27g| ) and ( p.27h| ) and considering the equations in the 



stationary limit. The terms ( p.27e ) and ( p.271|) are therefore in the following denoted as 
G-matrix-terms. 

In order to simplify the notation, the single particle basis is now chosen to diagonalize 
the unrenormalized mean-field Hamiltonian, i.e. 

haa' — ^aa'^a • (3.37) 

One can then rewrite the equation of motion for c in the limit described above as 

{lO — — efi) Cal3a'l3' — {uJ — e^/ — epi) Cal3a'l3' 

= E {^3X4\v\XiX2){Q {P2)\i\2a'f5' - <5ai A..','(P2).,A3A4} . 

A1A2A3A4 

(3.38) 

This equation is fulfilled for (r/ — > 0"^), if 
1 



ijj — e 



" - e/3 + A1A2A3A4 



E Qa/3A3A4^A3A4AiA2(P2o)AiA2Q'/3' (3.39) 
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and the conjugate equation hold, where at the right hand side the definition of the G- 
matrix has been inserted. Matrix-muhiphcation of ( p.39| ) with the matrix element of the 
two-body potential and eliminating c on the left hand side via c = p2 — P20 leads to 



(P20) 



20j7i72a'/3' 



7172 7172 



+ H {aP\v\-fl-f2) ^^^^^^^^^^ G'A3A4AiA2(P20)AiA2a'/3' 
7172A1A2A3A4 ^ ^71 ^72 + 

(3.40) 

for an arbitrary P20, i.e. we have 

Go^po.'p'= {aP\v\a'P')+ Yl MI^IAiAs) ^"^"^"^"^ . Gx,x,a'p' ■ (3.41) 

A1A2A3A4 uj - ea - + tr] 

Equation ( p.41| ) is a Bethe-Goldstone-equation for the G-matrix which guarantees the 
resummation of ring- and ladder- diagrams in vertical direction. 

3.5 Dynamical iteration in the horizontal direction 

The terms ( |3.27gD and (p.27h| ) each also contain a two-body correlation function coupled 
to a matrix element of the two-body potential, which leads to a dynamical resummation 
of diagrams. The topological structure of the terms, graphically represented in figs. |] 
and|^, in contrast to ( |3.27e| ) and ( |3.271| ) leads to a resummation in "horizontal" direction. 



where the factor Q-^ appears in the intermediate states and at a pair of an incoming and 
an outgoing index. 

In the single particle basis with paa' = ^aa'i^ai as for the terms accounting for the 
resummation in vertical direction, the summation over the intermediate indices breaks 
down and one remains with a factor of n'^ — Ua- In our notation each term containing one 
horizontal factor comprises a process together with its conjugate process. 



The direct term in ( |3.27g| ), with the indices in the matrix element of the two-body 
potential not interchanged, is shown in fig. |^ at the top. It leads to a dynamical resum- 
mation of ring-diagrams in horizontal direction. The exchange term in ( |3.27g| ) accounts 
for the dynamical resummation of ladder-diagrams in horizontal direction. 

In ( |3.27h| ) the direct term and the exchange term are equal; together they account for 
the resummation of ladder-diagrams in horizontal direction. 

All terms of the NCBCD approximation together lead to a dynamical resummation 
and mutual iteration of ring- and ladder-diagrams in vertical and horizontal direction on 
the two-body level, i.e. a parquet-like resummation pB| , |27[] which also renormalizes the 
one-body equation of motion. 



4 Numerical studies for pionic systems 
4.1 Comparison of partial resummations 

In the following we consider a pion gas as described in Appendix ^ confined to a box of 
sidelengths = Ly = = lOfm. In order to reduce the numerical effort we assume the 
pion gas to be in an isospin-symmetric configuration and thereby eliminate the internal 
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isospin degrees of freedom as described in Appendix p. As a single particle basis for our 
numerical simulations we choose standing waves, which vanish at the box boundaries. The 
calculations are performed using the 11 (in the noninteracting case) energetically lowest 
basis elements. 

Since the NCBCD equations only describe the propagation of the density matrices 
in time, we first have to solve the problem of finding proper initial conditions for the 
system at some given starting time. In the bosonic case there is no direct access to the 
manifold of stationary solutions of the NCBCD equations and, furthermore, it is e.g. not 
clear, what further constraints - besides the above mentioned sjTumetry properties - these 
solutions have to fulfill in order to remain physical. In |^ it was shown for the case of 
an analytically solvable model (Lipkin-model), that only requiring the density matrices 
to be positively definite and to satisfy the trace theorems for fermionic systems with good 
particle number leads to a sensible solution within the framework of the NQCD formalism. 
However, in our case the trace theorems do not hold. 

The most obvious starting point for the initial condition is given by a stationary self- 
consistent Hartree-Fock solution. In order to generate this solution for a given temperature 
and particle number, the one-body density matrix is occupied with a Bose-Einstein distri- 
bution with respect to the eigenstates of the mean-field Hamiltonian, where the chemical 
potential is adjusted in order to give the right particle number. 

In fig. H the Hartree-Fock energy levels and occupation numbers are shown for the 
example of systems with 60 and 90 pions for various temperatures, where the occupation 
numbers are given for only one isospin quantum number; the total occupation numbers 
are obtained by multiplication with a factor of 3. Due to the spatial symmetries of the 
system, some of the higher lying levels are degenerate. As one expects, the error due to 
the finite number of basis elements reduces in going to higher densities and lower thermal 
excitations. 

After providing a first approximation for the initial condition of the pion gas, we 
now investigate how the system evolves in time using the complete NCBCD equations or 
certain limiting cases for the partial resummation of diagrams on the two-body level. We 
aim at deciding, if for an approximately correct description of two-body quantities in the 
pionic many-body system one really needs to use the complete, numerically very involved 
NCBCD equations. 

The limiting case, in which only the resummation of ring- and ladder-diagrams in 
vertical direction is taken into account, is obtained from the complete equations ( |2.23|) 
and ( p.24|) by neglecting (|2.24d| ) and corresponds to time-dependent G-Matrix theory 
(TDGMT). It will be denoted by vertical approximation in the following. 

The limiting case, which only takes into account the resummation of ring- and ladder- 
diagrams in horizontal direction, is obtained from ( p.23| ) and ( |2.24| ) by neglecting ( p.24c| ) 



and will be denoted by horizontal approximation furtheron. In the fermionic case the 
approximation generated from the NQCD equations in an analogous way can be identified 
with a two-body RPA theory 0, |^. 

Using the isospin channel quantities defined in Appendix ^, we get the expression for 
the total energy of the system: 

E= {H) =3 5]Waa' (4.42a) 
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^ aa'f3l3' 

+ 1 E («'/3>s|«/?)3(p20)a/3a'/3' (4.42b) 

+1 {(^'(3>tu\a(3){3c^X,^,+6cX'p') 

^ aa'fifi' 

+^ E (3<*|„,^, + 6c^^„,^,) . (4.42c) 

The kinetic energy is given by (4.42a), the mean-field energy by ( [4.42b| ) and the correlation 
energy by (|442^ ); the total energy, which is conserved in time, is the sum of these three 
contributions. 

Even though in the investigation for fermionic systems in |^ the vertical as well as the 
horizontal approximation showed distinct deviations from the full NQCD approximation, 
the correlation energies still had about the same order of magnitude as in the complete 
correlation dynamics approach. For the pionic systems under investigation here this is no 
longer true. While at lower densities the 3 approximations do not differ significantly, at 
higher densities the horizontal approximation yields absolute correlation energies which 
are drastically higher than those extracted from the NCBCD calculations. 

In order to illustrate this behaviour, we show in the upper part of fig. |^ the correlation 
energies (averaged over time) for the 3 approximations as a function of the pion number 
at a temperature of 100 MeV. The system in each case was propagated up to the time 
t = 20 fm/c. One clearly sees the dramatic overestimation of correlation energy by the 
horizontal approximation with increasing density. In the lower part of fig. ^ we have 
left out the horizontal approximation in order to allow for a better comparison between 
NCBCD and the vertical approximation. The vertical approximation always generates a 
correlation energy, which is about a factor of 2 too small as compared to the NCBCD 
correlation energy. 

This demonstrates the importance of mixed diagrams consisting of contributions from 
the vertical and the horizontal approximation and of the interference of the different 
classes of diagrams. We are thus lead to the conclusion, that only the full NCBCD 
approach is able to yield adequate results for two-body quantities in the pionic system 
under consideration. Especially we infer, that a description of two-body dynamics by a 
Bethe-Goldstone approach or a time-dependent G-Matrix theory is not appropriate. For 
all further investigations we, therefore, use the complete NCBCD approach. 

4.2 Dynamical generation of a correlated state 

In this Sect, we will demonstrate, how one can dynamically generate an approximately 
stationary physical state for the NCBCD equations, which is fully correlated on the two- 
body level. 

A first attempt for the solution of this problem consists in initializing the Hartree- 
Fock solution and propagating the system in time until all correlations have built up and 
the system has reached a quasi- stationary state. However, in this state all density ma- 
trix elements will show more or less stable oscillations about an equilibrium value and 
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not necessarily approach a stationary state asymptotically, as we have checked numeri- 
cally. Furthermore, this approach has the additional problem, that the proper energy of 
the correlated state cannot be reached, since the NCBCD approximation conserves total 
energy. 

For the reasons mentioned above we introduce a different method by modifying the 
NCBCD equations of motion in the following way: The terms (2.23a) and (2.24a), which 
describe the free and mean-field propagation of p and c, remain unchanged. In all other 
terms, i.e. in (^^23b| ),( p^24b|) , (|^4|) and (|]24|), the matrix elements of the two-body 



potential {a(3\v\a' (3') are multiplied by a time- dependent dimensionless factor g{t). This 
alteration of the equations of motion, although at first sight appearing to be arbitrary, 
does not lead to an unphysical propagation of the system, since it is equivalent to an 
alteration of the Hamiltonian according to if — > H'{t), where 



H'{t) = J2 { taa' + (1 - 9{t)) E Px^xAaXMa'Xi) 



s 



A1A2 



+9i't) {c(f3\v\a'(3')al^alaa'af3'. (4.43) 

af3a'l3' 

li g(t) = 0, one recovers the TDHF approximation; if g{t) = 1 one recovers the full 
NCBCD approximation with respect to H. 

By initializing a stationary Hartree-Fock solution at time t = and propagating the 
system using the NCBCD equations induced by (|4.43|) up to the time t = tj, where 
g{0) = and g(tf) = 1, we continuously switch on the residual interaction on the two- 
body level. The total energy is not conserved during this process, thus enabling the 
system to energetically approach the correlated equilibrium configuration. 

Numerical experiments using different forms of g{t) have shown that a linear function 
requires the smallest amount of computation time in order to guarantee the convergence 
of the method. We, therefore, choose g(t) = Xt and hence tf = 1/A. 



In fig. |10] the change AE of the total energy reached at the end of the process is 
plotted versus tf for 30 pions in the volume of 10 fm^. The picture clearly shows, that 
there is a convergence as — 00 (adiabatic limit); i.e. if the parameter g{t) changes 
slowly, the response of the system is fast enough to generate an equilibrium configuration 
with respect to H'(t) at any given time t. The last statement implies, that for — > 00 
a convergence of the trajectory of the system in the overall configuration space has to 



occur. This is indeed the case, as can e.g. be seen from fig. where the correlation 
energy with respect to H during the process is depicted as a function of t/tf for various 
values of tf. Furthermore, fig. |rT] illustrates, that the oscillations caused by switching on 
the residual interaction too fast progressively vanish as one goes to higher values of tf, 
which clearly indicates that indeed in the limit — > 00 the system equilibrates at each 
time during the process. 

The change in total energy during the adiabatic process is not mainly caused by a 
buildup of correlation energy, but rather by a change in kinetic and mean-field energy by 
about the same order of magnitude, since the generation of two-body correlations causes 



a redistribution at the one-body level as well. This can be seen from fig. 0, where the 
changes in the single contributions to the total energy and the change of the total energy 
with respect to the Hartree-Fock solution are plotted as a function of tf. Furthermore, 
we note that for all densities and temperatures considered in this work, the total energy 
is lowered when the residual interaction is adiabatically switched on. 
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4.3 Two-body observables in coordinate space 

We now turn to the investigation of the two-body density matrix resulting from the 
adiabatic method described above. For this purpose the residual interaction in each case 
is switched on with t/ = 160 fm/c, and the values for p and c reached aXt = tf are used 
for the evaluation of the observables of interest. 

The two-body density matrix p2 in coordinate space is given by 

P2{rur2]r'ur'2) = E {p2)apa'p''^a{ri)^p{r2)^ljl,{r'i)^lj},{r'2) (4.44) 
and the uncorrelated part of p2, i-e. P20; and the two-body correlation function c by 

P2aifi,r2]r'i,r'2) = (p2o)a/3a'/3' ^a(n)^/3(^2)C'('"'l)#('^'2) (4.45) 

and 

C(ri,f2;r'i,r'2) = ^ C„^a'/3' V'a(n)^/3(^2)C'(^'l)^^'('"'2) , (4.46) 

with 

P2(ri, f2; r'l, r'2) = P2o(^i, ^2; r'l, r'2) + c(fi, f2; r'l, r'2) (4.47) 
and 

P2o(n, ^2; r'i,r'2) = p(ri; r'i)p(f2; r'2) p(fi; r'2)p(f2; r'l) . (4.48) 

The probability distribution for a simultaneous measurement of a particle at position ri 
and a particle at r2 is given by P2(^i, "^2; ri, r2). In fig. |1^ p(x, x; x, x) = / dydzp2{f, f; r, r) 
with r* = (x, I/, z) is plotted as a function of x for the temperatures 50 and 150 MeV 
at various pion densities (solid line). In addition, the same quantity is shown without 
the correlated part of p2, i.e. only for p2o (dashed line), and for the classical two-body 
distribution 

pf{ri, r2\ r'l, r'2) = p(ri; r"^i)p(f2; r'2) , (4.49) 

i.e. neglecting Bose statistics (dotted line). At the boundaries of the box all distributions 
vanish due to the boundary conditions imposed on the single particle basis elements. The 
inclusion of Bose statistics leads to an enhancement of the distribution obtained with 
P20 by a factor of 2 with respect to the distribution obtained with pg"**. Including the 
two-body correlation function leads to a further increase in the probability of finding two 
particles at the same place, which is due to the fact that the cr-interaction used is purely 
attractive. 

We now investigate the strength of the enhancement of p2{f,'r;f,f) relative to 
p{f] r) p{f] r) as a function of temperature and density. Naively one might expect, that 
with increasing density the system will be more correlated and therefore show a stronger 
enhancement. However, this is not true. In fig. |1^ one can already observe, that at 
equal temperature an increase in density does not necessarily imply an increase in the 
relative enhancement of the two-body probability. While at T=50 MeV the enhancement 
decreases with increasing density over the whole density range considered, at T=150 MeV 
there is a maximum in enhancement between 45 and 60 pions in (10 fm)^. 
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In order to illustrate this point more clearly the quantity 

pp J d^rp{f]r)p{f]r) 



(4.50) 



is plotted in fig. as a function of the pion density for various temperatures. For 
higher temperatures there are distinct maxima of the relative enhancement at certain 
densities. For all temperatures the relative enhancement then decreases again for higher 
densities. This behaviour is an indication of the nonperturbative nature of n-point cor- 
relation dynamics (the lowest order of which in our case is given by the Hartree-Fock 
approximation), since a denser and therefore more strongly interacting system does not 
imply a worse convergence of the cluster expansion. 

The next point to consider is the probability distribution for simultaneously finding 
two particles at different places in coordinate space. The question is, how this two-body 
density is altered by the inclusion of the two-body correlation function as a function of 
the relative distance of the two particles. 



In fig. 15 a two-dimensional cut through the six- dimensional quantities 



P2o(ri,f2;fi,f2) _ P2(ri,f2;ri,r2) 
p{fi]fi)p{f2;f2) ' p(ri;ri)p(f2;f2) 

at Hi = 1/2 = Zi = Z2 = 5 fm (i.e. in the center of the box of sidelength 10 fm) is depicted 
as a contour-plot as a function of Xi and X2 for two parameter pairs of density and 
temperature; in the left part of the figure only Bose-statistics are taken into account. The 
inclusion of the two-body correlation function (right part) leads to an increased gradient 
with respect to the relative distance, i.e. at smaller distances the two-body density is 
increased and at larger distances decreased. 



In fig. |I6| we show the same quantity as a function of x with xi = x and X2 = L — x, 
i.e. along a cut from the center to the lower right corner in the previous figure, for 50 
and 150 MeV temperature and various densities, where L = 10 fm is the sidelength of the 
box; the two-body densities evaluated with p2 (solid fines) and with p2o (dashed lines) are 
directly compared. 

Due to Bose-statistics, the curve for p2o always approaches 2 as x 0. Again the 
increased gradient with respect to the relative distance caused by the inclusion of the two- 
body correlation function can be seen, which is most pronounced for those parameters, 
that also show the strongest enhancement of probability at relative distance 0, i.e. at the 



maxima in fig. 14 



4.4 Relative strength of two- body correlations 

In the last Sect, we discovered, that for fixed temperature the relative enhancement of 
the probability for simultaneously detecting two particles at the same place in coordinate 
space exhibits a maximum at a certain density. This behaviour of the pionic many-body 
system will now be analyzed more closely. In this respect, we introduce the quantity 

_C_ ^ J2al3a'0' \ Ca/3a'l3'\ 

which will be denoted as "correlation strength" in the following. While in the two-body 
density in coordinate space only the isospin channel correlation functions c''*" contribute 
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(since we are looking at pairs of identical pions), all possible isospin configurations enter 
into the evaluation of the correlation strength ( [4.52[ ). Also all off-diagonal elements of 
the two-body density matrix are considered without having the possibility of mutually 
cancelling out each other. Although ( |4.52| ) is not invariant under unitary transformations 
of the single particle basis, it is still a useful measure of the relative importance of two- 
body correlations in the system. 

The curves in fig. |1^ were generated by adiabatically switching on the residual inter- 
action as described in Sect. ^]2| using tj = 160 fm/c and plotting the correlation strength 
reached aX t = tf versus density for various temperatures. Again there are distinct max- 
ima as in the case of the relative enhancement of the two-body density in coordinate 
space. If we compare fig. ^ to fig. we observe a very good agreement of the shape of 
the curves, implying that the observed result is not an artifact of the way we defined our 
measure for the relative importance of two-body correlations. 

As one approaches the maxima, the stationary Hartree-Fock solution used for the 
initialization at t = becomes increasingly useless. It is therefore tempting to assume, 
that the convergence of the cluster expansion in general might also become worse; an 
assumption which would have to be checked by explicitly taking into account higher 
order than the two-body correlation functions. However, within the framework of the 
NCBCD approximation it is not possible to verify this assumption. 

A possible interpretation of the observed maximum is that the system runs through a 
second order phase transition at a certain critical density (for fixed temperature), which 
corresponds to a critical chemical potential in the parameter space of the theory. In such 
a scenario one might expect the cluster expansion to break down as one approaches the 
critical region. Moving away from the critical region the mean-field (Hartree-Fock) solu- 
tion progressively becomes better and the convergence of the cluster expansion improves 



24|. 

A problem we want to point out in this context is that in a finite system, as considered 
in this work, one cannot draw any definite conclusions about the real order of a phase 
transition, since in such systems usually also first order phase transitions are smeared out 
and appear as second order phase transitions p9| . 

Finally, we explicitly demonstrate that not only a higher density, but also a stronger 
coupling does not automatically imply stronger correlations in the pionic system. For that 
purpose we show in fig. |T^ the correlation strength obtained by adiabatically switching on 
the residual interaction plotted versus density for 100 MeV temperature for the coupling 
given by ( |A.771 ) (solid line) and for a coupling twice as large (dashed line). The correlation 
strength for the higher value of the coupling constant is distinctly smaller than that for the 
lower one. This again clearly proves the nonperturbative nature of the cluster expansion. 



5 Summary 

In this work we have derived the NCBCD approximation for the nonperturbative, dy- 
namical description of particle number conserving bosonic many-body systems including 
two-body correlations. This approach leads to a dynamical resummation and mutual it- 
eration of ring- and ladder-diagrams in vertical and horizontal direction on the two-body 



level, i.e. a parquet-like resummation p6|, which also renormalizes the one-body 
equation of motion. 
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For our numerical simulations we assumed a pion gas confined to a box and a single 
particle basis of standing waves. The tttt interaction used was derived from a covariant 
Lagrangian in a three-dimensional reduction scheme and fitted to vrvr phase shifts. Within 
this model we showed, that a restriction of the equations to a resummation in horizontal 
direction leads to a dramatic overestimation of correlation energy and a restriction to a 
resummation in vertical direction leads to an underestimation of correlation energy by 
about a factor of 2. Thus an adequate description of dynamical two-body correlations 
in pionic systems requires the use of the complete NCBCD equations; e.g. a description 
within the TDGMT approximation is inappropriate. 

We, furthermore, showed a way to obtain a correlated, stationary and physical con- 
figuration by adiabatically switching on the residual interactions on the two-body level 
starting from a stationary Hartree-Fock solution. For this configuration we investigated 
two-body probability distributions in coordinate space and obtained an increased prob- 
ability for the simultaneous detection of two pions at the same place and a decreased 
probability for the simultaneous detection of two pions at large distances caused by the 
inclusion of the two-body correlation function. The two-body correlations in coordinate 
space as well as the correlation strength (defined in Sect. [4.4|) showed a distinct maximum 
in the relative strength of two-body correlations as compared to the disconnected parts 
of the two-body density matrix at a certain critical density at fixed temperature. This 
might be interpreted as a signature for a phase transition at a certain critical density. 
However, in the investigated temperature region of T < 200 MeV this maximum in each 
case was located at a density below 0.09/fm'^, which is much smaller than those estimated 
in ultrarelativistic nucleus-nucleus collisions |T0[| . 

Furthermore we found, that at higher densities the two-body correlations progressively 
decrease in relative importance, which implies the possibility of an adequate description 
of the system on the one-body level in such cases as ultrarelativistic heavy-ion collisions. 
We thus conclude that the conventional Hanbury Brown-Twiss (HBT) analysis of pion 
sources from vr — vr correlations in such reactions appears justified. 



Appendix 

A Specification of the model Hamiltonian 

In order to describe an interacting pion system we consider the Lagrangian 

£ = ^ {d,$) - ^m2$2 + i ^d,a) - ^MV^ - ga^' , (A.53) 

where $ = ($1, $2, '^'3) is the pion field of mass m and a is the scalar, isoscalar field of 
mass M mediating the interaction between the pions. 

From the Lagrangian (|A.53|) we get the coupled equations of motion for the classical 
fields or equivalently for the quantized field operators of the system: 

d''d^<l>i + rn^^i = -2ga^i (A.54) 

and 

d^'d^a + MV = -^$2 . (A.55) 



18 



In order to obtain a two-body interaction for the pions, we integrate out the cr-field. 
Formally inverting ( [A.55| ) we get 



a{k) 



It 



(A.56) 



With G„{k) 



k.k,-ivP+ie and G^{x) = f '^'^^G„{k) we have 



a{k) = gG^{k) J rfV e''=''"M$2(x') (i{x) = g / d^x'G^{x - x')<l>'\x 
Insertion of ( [A. 57] ) into ( |A.54| ) leads to 



d^d^ + m') <^i{x) = -2g'^i{x) / d^x'G^ix - x')^'{x') 



(A.57) 



(A.58) 



i.e. the cr-field has been eliminated from the equations of motion. 

We aim at a Hamiltonian, which yields ( [A.58| ) via the Heisenberg equations for $j and 
its canonically conjugate field momentum Ilj. 

This implies, that we have to carry out a three-dimensional reduction of the four- 
dimensional integration in ( A.58|) . As a first step we neglect the zeroth component of the 
four momentum fc^ in the denominator of (^^(/c), i.e. 



G^ik) ^ G^{k) 



- M2 



and 

GJx) 



d^k 



ikx 



-k^ - M2 



(A.59) 



(A.60) 



The intrinsic flaws of this naive instantaneous approximation will in part be cured later 
on. The easiest way to uncover these defects is to use a mode-expansion of the pion fields 
and to establish, which terms are responsible for s-, t- or u-channel scattering. 
With (|A39|) and ( [OOD we get 



d^x'G^{x-x')$\x' 
d^k 



d'^k 
2^ 



-^^ix" 



- M2 



d^x' 







d^x'G^{x-x')$\x',x^) . 



6{x^ - X' 



(A.61) 



The equation of motion from Hamiltonian dynamics in the instantaneous approxima- 
tion (|A.61|) thus reads: 



[d^d^ + m^) = -2/$i(f) J d^x'G^ix- x')$^{x) , 
where all field operators are taken at the same time t. This is fulfilled for 
H = ^J d^x |n^(x) + (V$(f))^ + m2$2(a 



(A.62) 



+^g'' I d:^xd^x'^\x)G„{x - x )<l>'\x 



(A.63) 
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since (|A.63|) together with ( p.6|) and the canonical equal-time commutation relations 

= [n,(f,t),n,(f' = , 

[<^i{x,t),Uj{x',t)] = i6ij6^^\x - x) (A.64) 
generate the equations of motion 

dtM^) = Ui{x) (A.65) 
and 

(A.66) 



dtUi{x) = V^<l>i(f) - m2$i(f) - 2c/2$i(f) J d^x'G^ix- x')<l>\x) , 

equivalent to ( |A.62| ). 

The mode expansions for the pion field operators now read: 



and 

ni(f) 

with 



Akx 



(27r)2 ^2uj{k) 



a{k, i) + a){—k, i) 



(A.67) 



d^k 



Akx . 



u{k) 



a{k, i) — a\—k, i) 



k^ 



uj{k) = 
and 

a{k,i),a\k',j)] = 6ij6^^\k-k') 



(A.68) 

(A.69) 
(A.70) 



where i denotes the isospin of the corresponding operator with regard to the cartesian 
representation of SU(2). Inserting the mode expansions ( |A.67| ) and ( |A.68| ) into ( |A.63| ) 
leads to 



H = Y,jd^k u{k)a\k, T)a{k, t) 



+ J2 d^kid^k2d^k3d^k4- 
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(A.71) 

where now r denotes the isospin quantum number. 

The Hamiltonian ( |A.71| ) still contains terms (underlined) not commuting with the 
particle number operator. In view of Sect. ^ we neglect such terms in our approach, 
i.e. we only consider elastic scattering of pions (also including particle-number conserving 
s-channel scattering). This is equivalent to the assumption, that the - particle number 
violating - operator products generated by the underlined terms in ([A.71|) (after normal 



ordering) in the equations of motion for particle number conserving operator products 
have vanishing or at least negligible expectation values, in agreement with assumption 
(|A|) of Sect. |. 

Normal ordering the remaining terms in ( |A.71| ) gives: 



H = Y, j £k uj{k)a\k, T)a{k, t) 

2 

x5^^\ki + k2 + h + h)Ga{ki + h) 

X |a^(— ^3, r')a^(— A;4, T')a{ki, T)a(k2, t) + a\—k2, T)a^(— A;4, T')a{ki, T)a(k3, r') 
+a\—k2, T)a^(— A;3, T')a{ki, T)a(k4, r') + a\—ki, r)a^(— A;4, T')a(k2, T)a{k3, r') 
+a\—ki,T)a\—k3, T')a{k2, T)a{ki, r') + a){—ki, T)a){—k2-, 'T)a{kd,i T')a{k4, r')| 

= ^ / dPk uj{k)o){k,T)a{k,T) 

+ Y.j d^k^d^k2d^k3d''k^-^ {ujihMk2)ujih)ioih)y~' S^^\h + k2-h- h) 

X {2(5^(^1 + ^2)0^(^1, 'r)a^{k2, T)a{ks, T')a{ki, r') (A. 72a) 

+4G,(A^i-4)a^(A?i,r)a^(4,r')a(A?3,r)a(fc4,r')} • (A.72b) 



As advertised before, one can now extract the terms corresponding to the different 
scattering channels; i.e. (A. 72a) is responsible for the s-channel scattering of pions, since 
the intermediate cr-propagator depends on the sum of momenta in the in- or outgoing 
channel and the (cartesian) isospin quantum numbers of the two pions in each of these 
channels have to be equal. 

Term ( |A.72b|) is responsible for t- and u-channel scattering of pions, since the inter- 
mediate (T-propagator depends on the momentum transfer from one of the incoming to 
one of the outgoing particles, while these two particles have to carry the same isospin 
quantum number as well as the remaining two pions. 

We are now in the position to improve the three-dimensional reduction scheme. In 
the center of momentum-frame of two colliding pions the zeroth component of the four 
momentum in the a-propagator vanishes for t- or u-channel scattering, while for s-channel 
scattering the three momentum in the cr-propagator vanishes. This implies that the naive 
instantaneous approximation ([A.59| ) is inappropriate for the s-channel term of the two- 
body potential. 
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The situation is now cured by simply placing the zeroth components of the pion 
four-momenta on-shell. The ambiguity, which of the - up to now equivalent - three mo- 
mentum combinations in the a-propagator have to be used, we resolve by averaging over 
both possibilities. This method is in close analogy to the well-known three-dimensional 
reduction schemes of Gross, Blankenbecler-Sugar and Thompson except for the dif- 
ference, that with these methods one carries out the three-dimensional reduction of T- 
or G-matrix equations, where in addition to the external three momenta the invariant 
mass of the colliding system is fixed. In the latter case one has different possibilities 
of placing the particles off-shell, such that their four momenta give the correct ^/s. In 
a dynamical many-body calculation as in this work, the y/s for an elementary two-body 
scattering process can only be reconstructed from the three-momenta involved, so that 
the schemes cited above are not directly applicable. 

Since point-like interaction concepts are inappropriate for hadron scattering we in- 
troduce formfactors for the TTTrcr-vertices containing momentum cutoffs to regularize the 
theory. We use a formfactor of the form 

where A is the cutoff parameter and k is the four momentum of the cr-particle coupling 
to the vertex. The formfactor (|A.73|) has the distinct advantage of leading to a uniform 
parameterization for spacelike and timelike four momenta without running into a pole in 
one of the two domains. 

The Hamiltonian we use in our explicit computations then reads: 

H = J2 fd^k uj{k)a\k, T)a{k, r) 

+^ E / d^hd^k2d^hd^h 5^-^\ki + h - h - h) 

tt' 

X [{hk2\vs\ksk4)a'^ {ki, T)a\k2-, T)a{ks, T')a{k4, r') 

+ {kik2\vt,u\k3k4)a\ki, T)a\k2, r')a{k-i, T)a{ki, r')| (A.74) 
with 

-.2 



{kik2\vs\hk4) = ^27r)3 (^(^i)^(^2)t^(4)t^(fe 
{F{iujih)+u;{k2)?-{h + k2)')y 



2 



[F{{u;{h)+uj{k4)f-{h + hf)} 

and 

f 

{f {(ujih) - ujih)? - {h - h?)} 



(A.75) 



{kihlvt^ulhh) = ^7^-^ {uj{ki)uj{k2)uj{k3)uj{k 



4. 

2 



X 



(uJ{k^) - Lu{k,))^ - (k^ - ksY - 
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(iuih) ~ iuihW - ih - - ^ 

The a can be associated with the resonances in the {J^'" = 0++, I'^ = 0+)-channel 
(cr-channel) in Trvr-scattering, where the theory described can easily be generahzed to a 
couphng of the pion field to more than a single scalar and isoscalar field; we simply obtain 
the matrix elements of the two-body potential as a sum over those of the form ([A.75| ) 



( |A.76 ) for the individual a-particles (jj of masses Mj with couplings Qi and cutoffs Aj. The 



parameters used in this work are taken from pl[]p|, where additional particles besides the 
scalar and isoscalar ones are used to fit scattering data, but the couplings for the latter 
are directly deduced from the widths of the resonances in the cr-channel by means of a 
Breit-Wigner formula. Using only cr-mesons - in order to keep the numerical task within 
a reasonable range - should at least guarantee a realistic description of the short-range 
attractive part of the vrvr-interaction. The parameters taken from pT| are: 



Ml = 980 MeV, Ai = 1200 MeV, gi = 595 MeV 
Ma = 1300 MeV, A2 = 1200 MeV, c/2 = 1854 MeV 

m =m^ = 140 MeV . (A. 77) 

In view of the complexity of the NCBCD equations we do not claim to have a fully 
realistic model of a pion gas, which would require the inclusion of p-mesons and kaons in 



a coupled channel-calculation |32|-|38| in order to fit the experimentally measured phase 
shifts for the higher partial waves [^, The interaction should work at moderate 
relative momenta. 



B Equations of motion for the isosymmetric 
pion system 

Since for a numerical solution of the NCBCD equations the required computation time 
increases like and the required memory space increases like N'^, where N denotes 
the number of single particle basis states, it is furthermore desirable to eliminate all 
internal degrees of freedom from the equations and thereby reduce iV by a factor of 3. We 
therefore assume, that the system is in an isospin-symmetric state and that all density 
matrices are diagonal in isospin space, which in the spherical representation of isospin 
SU(2) implies, that the sum of the isospin-z-components in the incoming and outgoing 
channels have to be equal. As one easily verifies, isosymmetry and isodiagonality are 
conserved dynamically; therefore these assumptions do not lead to an inconsistent theory. 

Now we use a = {a, To), where a denotes the quantum numbers of the space-time 
degrees of freedom in an arbitrary orthonormal single particle basis and r denotes the 
isospin quantum number in the cartesian representation. The Hamiltonian of the pionic 
system can then be written as 

H = '^Ua'a]^aa' + ^ {a(3\v\a'(3')al^a^^aa'a^, (B.78) 

^We use a slightly different form of the formfactor, which however agrees reasonably well with their 
form in the kinematical regime under investigation. 
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with 



taa' — taa'^TaT^ (B.79) 

and 

{aMo^'~P') = {a(5\vs\a'(5')5r^r^5r^,r^, + {a(5\vt,u\a' (5')8r^rjr,r,, ■ (B.80) 

Isosymmetry implies, that in the cartesian representation of SU(2) all density matrices 
have to be invariant under an overall permutation of isospin quantum numbers. 
For the one-body density matrices we assume 

Paa' — Paa'^TaT^, ■ (B.81) 

For the two-body correlation function c-^^,^, there are 4 possible distinct isospin con- 
figurations: 

1- Ta = Tf3 = Ta' = Tpt 

2. Ta = T/s , Tq,/ = rg/ , Tq, ^ Ta' 

3. Ta = To,' , Tp = Tp> , Tq, 7^ 

4. Ta = Tpi , T0 = Ta' , Ta T^ . (B.82) 

All other configurations, which do not contain two pairs of equal isospin indices, violate the 
diagonality of the two-body density matrix in isospin space; the corresponding correlation 
functions are thus assumed to vanish. 

Configuration 1. corresponds to the situation, where s-, t- and u-channel scattering 
contribute to the Born amplitude for the elementary process; configuration 2. corresponds 
to the situation, where only s-channel scattering contributes and configurations 3. and 4. 
reflect the situation, where only t-channel scattering contributes. 

We can therefore decompose the two-body correlation function in the following way: 

'^apa'P' ~ ^TcTp^T^iTpi^TcT^iCajSa' P' + ^TcTb^t^iTi^X^ ~ ^TaT^,)Cafia' jl' 

+^TaT„/^r^T^/(l — ^raT^)Ca^a'/3' + ^TaT^/^r^r^, (1 " ^TcTi3)c^al}a' P' (B.83) 

with the isospin channel correlation functions c''*", c*, c*^, c*^. 
The one-body density matrices fulfill 

Paa' = {pa' a)* , (B.84) 

while the isospin channel correlation functions fulfill 

stu,s stu,s _ stu,s _ ( stu,s \* ^'"R Sc;'\ 

^aPa'/3' — '-paa'p' — '^ap/i'a' — \'-a'/3'a/3) {D.QO) 

and 

'^a/3a'f3' = '^faa'/S' = <^a/3/3'a' = i^a'fS'a/s) 5 C^^/3q'/3' = i'^a'fS'als) ■ (B.86) 

We note in passing, that for our purpose the cartesian representation of SU(2) is 
technically superior to the spherical one, because in the spherical representation we have 
to consider 7 instead of 4 different isospin channel correlation functions. 

We can now explicitly carry out all summations over isospin indices appearing in the 
NCBCD equations, which leads to a coupled system of a single one-body equation of 
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motion and four two-body equations of motion, where all considered quantities now only 
depend on space-time degrees of freedom. 
With the abbreviations 

Uaa' = Yl ^AiA2 (3(Q;A2|vt,„|Q;'Ai) + (aAalvt^^l AiQ;') + {aX2\vs\a'Xi)g) , (B.87) 
A1A2 

haa' = taa' + Uaa' , (B.88) 







(B.89) 
(B.90) 



we obtain the NCBCD equations for the isosymmetric, isodiagonal pionic system: 

'idtPaa' = XI i^aXPXa' — hxa' Pax) 
X 

+ Y: {(«Ai|^,„|A2A3) (Ca3,,,, + 24i,3„,,J + (aAil^.lA^As) (^3,,,, + 2c^^,3,,,^_ 
A1A2A3 

-{X,X,Ma'X,) {cXm^^ + ^c'Xmx.) - (AiA2k>'A3) (C,, + 2<,3,,,J} (B.91) 



and 



idt 



~ ^stu 
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^stu 

'^0/3x3x4 
^0^X3X4 

^a/SXsXi 



- {X3X4\Vs\XlX2)QxM, 
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r stu I 9 tl 



„stu I 9„tl 
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3 Figure Captions 

fig.l : The Hartree-terms in the t-, u- (l.h.s.) and s-channel (r.h.s.); dashed hori- 
zontal hnes are spacehke, dashed vertical lines are timelike cr-fields, solid lines 
represent pions. 

fig. 2 : The Fock-terms in the t-,u- (l.h.s.) and s-channel (r.h.s.). 

fig.3: The term coupling to the two-body correlation function; t-, u-channel (l.h.s.), 
s-channel (r.h.s.). 

fig.4: Contributions to the Born-term. 

fig.5: Resummation terms in vertical direction with in the incoming channel; 
t-,u-channel term (l.h.s.), s-channel term (r.h.s.). 

fig.6: Resummation terms in horizontal direction in the t-, u-channel; direct term 
(upper part), exchange term (lower part). Note that connects a and A3, 
a' and Ai. 

fig.7: Resummation terms in horizontal direction in the s-channel; direct and ex- 
change term are identical. 

fig.8: Single particle energy levels and occupation numbers of the Hartree-Fock ini- 
tialization for 60 (upper part) and 90 (lower part) pions in a volume of (10 fm)^ 
at various temperatures. The symbols denote the position of the discrete en- 
ergy levels. 

fig.9: Time-averaged correlation energy as a function of the pion number at 100 
MeV temperature for NCBCD, vertical and horizontal approximation (upper 
part); magnification (lower part). 

fig. 10: Change in total energy after switching on the residual interaction as a func- 
tion oi tf. 

fig.ll : Correlation energy with respect to if as function oit/tf for different values 

fig. 12 : Change of the single contributions to the total energy as function oi tf; A 
kinetic energy, B mean-field-energy, C total energy, D correlation energy. 

fig. 13: / dydzp2{f,f;r,f) (solid line), / dydzp2o{r,r;r,r) (dashed line) and 

/ dydzp{f;r)p{f;r) (dotted fine) as a function of x at 50 MeV (upper part) 
and 150 MeV (lower part) temperature for various densities. 

fig. 14: / d^rp2{r, r; r, r)/ / d^rp{r; f)p{f; r) for various temperatures as a function of 
density. 

fig. 15: Cut through the two-body density divided by the product of one-body den- 
sities at = ^2 = -^1 = -^2 = 5 fm at 50 MeV temperature and 30 pions 
in (lOfm)^ (upper part) and at 200 MeV and 90 pions (lower part); P20/PP 
(l.h.s.), P2/PP (r.h.s.). 



29 



fig.16: One-dimensional cut through the two-body density divided by the product 
of one-body densities at 50 MeV (upper part) and 150 MeV (lower part) 
temperature and for various densities; the relative distance is parameterized 
by X. Solid line P2/ PP, dashed line pi^j pp. 

fig. 17: Relative correlation strength after switching on the residual interaction with 
tf = 160fm/c. 

fig. 18: Correlation strength at 100 MeV temperature; solid line with normal cou- 
pling, dashed line with doubled coupling. 
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